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ABSTRACT 

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical 
observations. We have recently determined accurate values for its mass and distance, and for the orbital incli- 
nation angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical 
model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal 
continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the 
black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains 
a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3a). For a less probable (synchronous) 
dynamical model, we find a* > 0.92 (3a). In our analysis, we include the uncertainties in black hole mass, 
orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux 
via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced 
by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole 
and the disk's low luminosity. 

Subject headings: accretion, accretion disks - black hole physics - stars: individual (Cygnus X-1) - X-rays: 
binaries 



1. INTRODUCTION 

In our two companion papers (Reid et al. 201 1 ; Orosz et al. 
201 1), we report accurate measurements of the distance, black 
hole mass and orbital inclination angle for the black-hole bi- 
nary system Cygnus X-1. Herein, we use these results to de- 
termine the spin of the black hole primary by fitting the ther- 
mal component of its X-ray spectrum to our implementation 
of the Novikov-Thorne model^ of a thin accretion disk (Li 
et al. 2005). 

Cygnus X-1 was discovered at the dawn of X-ray astron- 
omy (Bowyer et al. 1965) and is one of the brightest and most 
persistent celestial X-ray sources. Its compact primary was 
the first object to be established as a black hole via dynami- 
cal observations (Webster & Murdin 1972; Bolton 1972). For 
several decades, the source has been extensively observed at 
radio, optical, ultraviolet and X-ray wavelengths. 

Cygnus X-1 is typically in a hard spectral state, but occa- 
sionally it switches to a soft state, which may persist for up 
to a year (see Figure 1). It is only in this soft state, when 
the disk spectrum is prominent, that one can measure the spin 
using the continuum- fitting method we employ. In this soft 
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State, which corresponds to the steep power-law (SPL) state^, 
a strong Compton component is always present. Although 
Cygnus X-1 has been observed on thousands of occasions, it 
has never been observed to reach the thermal dominant (TD) 
state, the state that is most favorable for the measurement of 
spin via the continuum- fitting method (Steiner et al. 2009a). 

Following the pioneering work of Zhang et al. (1997), we 
measure black hole spin by estimating the inner radius of the 
accretion disk Ri^. In this approach to measuring spin, one 
identifies Rin with the radius of the innermost stable circular 
orbit i^isco, which is predicted by general relativity. Risco 
is a monotonic function of the dimensionless spin parameter 
a*, decreasing from 6 GM/c^ to 1 GM/c^ as spin increases 
from a* = Otoa* = 1 (Bardeen et al. 1972)^. This rela- 
tionship between a* and i^isco is the foundation for measur- 
ing spin by either the continuum- fitting method (Zhang et al. 
1997) or by the Fe Ka method (Fabian et al. 1989; Reynolds 
& Nowak 2003). 

The identification of Rin with Risco is strongly supported 
by the abundant empirical evidence that the inner radius of 
the disk in the soft state of a black hole binary does not ap- 
pear to change even as the temperature and luminosity change. 
LMC X-3 is a prime example; its inner-disk radius has been 
shown to be stable to several percent over a period of 26 years 
(Done et al. 2007; Steiner et al. 2010). Strong theoretical 
support for identifying Ri^ with Risco is provided by mag- 
netohydrodynamic simulations of thin accretion disks, which 
show the disk emission falling off rapidly inside the ISCO 
(Shafee et al. 2008; Reynolds & Fabian 2008; Penna et al. 
2010; Kulkarni et al. 2011; Noble et al. 2010). 



^ Throughout, we use the X-ray states defined by Remillard & McClin- 
tock (2006): hard; thermal dominant (TD); steep power-law (SPL); and in- 
termediate states. In the alternative state classification scheme of Roman & 
Belloni (2005), the only states reached by Cygnus X-1 are the Low/Hard, 
Hard-Intermediate and Soft-Intermediate, which correspond respectively to 
our hard, intermediate (i.e., hard:SPL) and SPL states. 

^ a* = cJ/GM'^ with |a*| < 1, where M and J are respectively the 
black hole mass and angular momentum. 
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In our early work on measuring spin, we relied solely on 
TD-state data (e.g., Shafee et al. 2006; McClintock et al. 
2006; Liu et al. 2008). More recently, using our empirical 
model of Comptonization SIMPL (Steiner et al. 2009b), we 
have shown that one can obtain values of the inner disk radius 
using SPL data that are consistent (within ~5%) with those 
obtained using TD data if the scattering fraction fsc ^ 25% 
(Steiner et al. 2009a). Consequently, we are now able to rou- 
tinely and reliably obtain spin estimates for sources in the SPL 
state (Gou et al. 2009; Steiner et al. 2011) if the scattering 
fraction is not too extreme. This development has paved the 
way for measuring the spin of Cygnus X-1, which has never 
been observed in the TD state. 

For the continuum-fitting method to succeed, it is essen- 
tial to have an accurate value for the black hole mass M in 
order to express the radius appropriately as the dimension- 
less quantity Riscoc^ /GM. For the method to succeed, it is 
furthermore essential to know the disk luminosity, and there- 
fore to have accurate values for the distance D and the disk 
inclination angle i, which we infer by assuming that the or- 
bital angular momentum vector and black-hole spin axis are 
aligned (see Section 7.4). We measured these three critical pa- 
rameters in a two-step process: Firstly, using the Very Long 
Baseline Array (VLB A), we determined a model-independent 
distance via trigonometric parallax that is accurate to 6%: 
D = l.Seto n kpc (Reid et al. 2011). Secondly, using this 
accurate distance estimate to constrain the radius of the com- 
panion star, we then modeled an extensive collection of opti- 
cal data. For our favored asynchronous model. Model V, we 
find: M = 14.8 ± 1.0 and i = 27A± 0.8 deg (Orosz 
et al. 2011). With the values of these three input parameters 
in hand, and using three soft- state X-ray spectra, we deter- 
mined the spin of Cygnus X-1; our final value is presented in 
Section 6. 

Spin measurements have previously been obtained using 
the continuum-fitting method for eight black holes (Mc- 
Clintock et al. 2011). Three have low to moderate spins, 
a* < 0.5; four have high spins, a* ~ 0.7 — 0.9; and one 
of them, the archetypal microquasar GRS 1915-1-105 (Mirabel 
& Rodriguez 1994; Fender & Belloni 2004), is a near-extreme 
Kerr hole with a* > 0.98 (McClintock et al. 2006; Blum et al. 
2009). 

The paper is organized as follows. In Section 2, we dis- 
cuss our observations and the selection and reduction of the 
X-ray spectral data. The analysis of these data and our re- 
sults are presented respectively in Sections 3 and 4 for our 
adopted model, and the robustness of these results is discussed 
in Section 5. (Additional analysis work and results for several 
preliminary models are presented in the appendices.) In Sec- 
tion 6, we first determine the error in the spin parameter due 
to the combined uncertainties in D, M, i and the absolute flux 
calibration, and we then present our final value of a* and con- 
fidence limits. Seven distinct topics are addressed in Section 
7, and in Section 8 we offer our conclusions. 

2. DATA SELECTION, OBSERVATIONS AND DATA REDUCTION 

There exist very few spectra of Cygnus X-1 that are suit- 
able for the measurement of the black hole's spin for reasons 
we now discuss. A typical soft-state (and SPL) spectrum is 
comprised of three principal elements: a thermal component, 

fsc is the fraction of the thermal seed photons that are scattered into 
the power-law tail, and it is the normalization parameter of our XSPEC model 
SIMPL/SIMPLR (see Section 3.2). 



a power-law component, and a reflected component, which 
includes the Fe Ka emission line. Three such spectra are an- 
alyzed in detail and illustrated in Section 3. It is apparent 
from an inspection of these spectra that the spectral cover- 
age must extend to at least 30 keV in order to constrain the 
strong power-law and reflection components. At the same 
time, because the temperature of the thermal component is 
consistently low (kT ^0.5 keV; Appendix A), one also re- 
quires coverage down to 1 keV in order to constrain this 
crucial thermal component, which is partially absorbed at low 
energies by intervening gas. The rarity of spectra that meet 
these requirements became clear to us following our exhaus- 
tive search of the thousands of spectra of Cygnus X-1 that are 
contained in the HEASARC data archive. To our surprise, we 
found only a single spectrum, SPl (Table 1), that is suitable 
for the measurement of black hole spin^^ It was obtained 
on 30 May 1996 in an observation made simultaneously us- 
ing the Advanced Satellite for Cosmology and Astrophysics 
(ASCA) and the Rossi X-ray Timing Explorer (RXTE). 

There is a second, important reason for the paucity of suit- 
able spectra, namely that Cygnus X-1 is seldom in the re- 
quired disk-dominated state. This fact is illustrated in Fig- 
ure 1, which summarizes the behavior of Cygnus X-1 since 
1996 as observed using the RXTE AU-Sky Monitor (ASM). 
As indicated in the figure, we select only those data for 
which the spectral hardness SH < 0.7, which occurs < 
10% of the time. Fortunately, Cygnus X-1 entered its soft 
state in mid-2010, and we obtained two additional broad- 
band spectra on 22 July and 24 July 2010 by making simul- 
taneous observations using the Chandra X-ray Observatory 
and RXTE (Table 1). The times of these two observations 
and the ASCA/RXTE observation are indicated by arrows in 
Figure 1, and the corresponding ASM measures of intensity 
and spectral hardness are plotted as red stars. Detailed in- 
formation on these three observations is summarized in Ta- 
ble 1. Three spectra were derived from these observations: 
the 1996 archival ASCA/RXTE spectrum SPl, and the two 
Chandra/RXTE spectra, SP2 and SP3, which were obtained 
in 2010. 

We first focus on spectrum SPl, which has been described 
in detail by others (Dotani et al. 1997; Cui et al. 1998; 
Gierlinski et al. 1999). We include in our analysis all of the 
data collected by both ASCA and RXTE, even though the ob- 
servations were strictly simultaneous for only about 1000 s. 
This is a reasonable approach because both missions show 
that the source intensity was stable during the entire observing 
period (see Figure 4 in Gierlinski et al. 1999). For A^CA, we 
consider only data collected by the second Gas Imaging Spec- 
trometer (GIS) detector, GIS2 (the GIS3 data were excluded 
because of an unexplained residual feature in the spectrum 
near 1 keV). We disregard the Solid-state Imaging Spectrom- 
eter (SIS) data because of data-rate limitations. For RXTE, we 
use all of the data collected by all five Proportional Counter 
Array (PCA) detectors. For the main analysis, we disregard 
the RXTE High Energy X-ray Timing Experiment (HEXTE) 
data because the useful bandwidth of the PCA, which extends 
to 45 keV, already provides more than adequate energy cover- 
age. We demonstrate this fact in Section 5.2 where we show 

We decided not to use the Swift XRT/BAT spectra, which provide 
broadband coverage, because (1) the XRT timing spectra (e.g., Swift ID 
00031651005) show strong residual features in the energy range 1-3 keV; (2) 
the Low-Rate Photodiode spectrum (Swift ID 00101469000) is poorly caH- 
brated (Romano et al. 2005); and (3) Swift provides no coverage at all from 
10-15 keV. 
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that the inclusion of HEXTE data has no significant affect on 
our results. 

The GIS spectra were extracted following the standard pro- 
cedures described in The ASCA Data Reduction Guide^^. The 
GIS2 spectrum was fitted over the energy range 0.7-8.0 keV. 
We did not correct the effective area of the detector because 
the instrument team had already done this calibration to an 
accuracy of 3% using the standard Crab spectrum (Toor & 
Seward 1974; Makishima et al. 1996; Steiner et al. 2010). 

We reduced the RXTE PCA data following the same pro- 
cedures described in McClintock et al. (2006). Data reduc- 
tion was performed with standard tools from the HE AS OFT 
package provided by NASA. The critical steps in determining 
the PCA background and making the response files for spec- 
tral analysis utilized HEASOFT version 6.10. X-ray spectra 
were extracted from the Standard 2 telemetry mode, which 
provides coverage of the full PCA bandpass every 16 s. Data 
from all Xe gas layers of PCU-2 were combined to make each 
spectrum. The background spectrum was determined with the 
"bright source" model. Redistribution matrix files and an- 
cillary response files were freshly generated individually for 
each PCU layer and then combined into a single response file 
using the tool pcarsp. We used the PCA response matrices^-^ 
vl 1 .7 released on 2009 August 17, which allowed us to obtain 
reliable fits over the energy range 2.55-45.0 keV. 

We corrected for the effective area of the PCA using the 
spectrum of the Crab Nebula as a standard source and us- 
ing our recently-adopted and improved method, which is de- 
scribed in Steiner et al. (2010). Specifically, we compared 
the Crab spectrum of Toor & Seward (1974) (T = 2.1 and 
N = 9.7 photons keV~^) to parameters obtained by an- 
alyzing proximate archival observations of the Crab. In this 
way, we determined a pair of correction factors for spectrum 
SPl: a normalization correction factor Cts = 1.123 (the ra- 
tio of the observed normalization to that of Toor & Seward) 
and a correction to the slope of the power law, ATts = 0.023 
(the difference between the observed value of F and that of 
Toor & Seward). These corrections were applied in all of 
our analysis work via a custom XSPEC multiplicative model 
CRAB COR. We also corrected the detector count rates for 
dead time by the factor 1.048. 

We turn now to the recent Chandra/RXTE observations (Ta- 
ble 1; spectra SP2 and SP3). The observations performed by 
the two spacecraft were mostly simultaneous, overlapping for 
about 4 ks, while the total duration of each Chandra obser- 
vation was ^ 6 ks. Because the source intensity for both 
observations was constant to within 10%, we included in our 
analysis all of the Chandra data, plus the strictly simultaneous 
RXTE FCA data. 

Our Chandra observations were made using the High- 
Energy Transmission Grating (HETG) and the Advanced 
Camera for Imaging and Spectroscopy (ACIS) (Canizares 
et al. 2005; Garmire et al. 2003). The data-rate limitation of 
the detectors makes it challenging to observe a bright and vari- 
able source like Cygnus X-1. The principal problem is "pile- 
up," i.e., the registering of two or more photons in the same 
or adjacent pixels within a single frame time. Given the un- 
certainties, and as a hedge, we performed the pair of Chandra 
observations using two different instrumental configurations, 



Timed Exposure (TE) mode and Continuous Clocking (CC) 
faint mode, which have contrasting virtues and limitations. 

In reducing the Chandra TE-mode data, we followed the 
method described by Smith et al. (2002). In order to avoid 
saturating the telemetry, only the data for the High Energy 
Grating (HEG; -1 order) and Medium Energy Grating (MEG; 

order) components of the HETG were recorded. The data 
for the readout streak on the same side of the HEG and MEG 
spectra were also recorded. We extracted the "streak" and 
background spectra following the recommended procedures 
Although the net exposure time for the HETG spectrum is 
2.15 ksec, the effective exposure time for the streak spectrum 
is only about 8 sec^^. For this spectrum, we estimate that less 
than 3% of the events are affected by pile-up, and we therefore 
use the full 0.5-10 keV bandwidth. For the dispersed grating 
spectrum, we only included data for which < 5% of the events 
are piled up; for the HEG and MEG respectively these data are 
in the energy ranges 0.7-0.9 keV and 7.0-10.0 keV. 

In the CC mode, the frame time is reduced to about 3 ms 
(compared to 1.3 s for the TE mode) by continuously trans- 
ferring the data from the image array to the frame- store ar- 
ray. While this essentially eliminates pile-up, the details of 
the spatial distribution in the columns are then lost due to the 
collapse of the 2D image into a ID image. Again, only the 
orders of the HEG and MEG spectra mentioned above were 
recorded. The spectra were extracted following the standard 
procedures'^. We fitted these data over the full energy range 
0.5-10.0 keV, except for the 1 .3-2.0 keV chip gap in the MEG 
spectrum. A downside of the CC mode is that its calibration 
is less certain than that of the TE mode. 

For RXTE, lacking data from all five Proportional Counter 
Units (PCUs), we elected to use only the data from what his- 
torically has been the best-calibrated detector, PCU2. (One 
obtains essentially the same results for the Chandra/RXTE 
spectra using any combination of the available PCUs; like- 
wise, for the ASCA/RXTE spectrum, it is unimportant whether 
one uses all PCUs, as we did, or PCU2 alone.) As described 
above, we corrected the effective area of the detector. For both 
observations, the normalization correction is Cts = 1.073, 
and the power-law slope correction is AFxs = 0.029. The 
dead time corrections for SP2 and SP3 are respectively 1.052 
and 1.044. 

Finally, we included systematic errors in the count rates in 
each PHA channel to account for uncertainties in the instru- 
mental responses: 1% for GIS2 and 0.5% for all the PCUs. 
All the Chandra data were binned to achieve a minimum num- 
ber of counts per channel of 200; no systematic error was in- 
cluded because the statistical error is large. 

3. DATA ANALYSIS 

A typical spectrum of Cygnus X-1 is comprised of three 
principal elements: a thermal component, a power-law com- 
ponent, and a reflected component that includes the Fe Ka 
emission line. The structures in the X-ray source that give rise 
to these components, namely the accretion disk and its corona, 
are illustrated in Figure 2. The spectral components them- 
selves, in relation to the total observed spectrum, are shown 
plotted in Figure 3, which illustrates the results of the rela- 
tivistic analysis described in this Section. 



1^ The manual is available at http://heasarc.gsfc.nasa.gov/docs/asca/abc/abc.html http://cxc.harvard.edu/ciao/threads/streakextract/ 

For a description of the latest response files, see The procedure for calculating this exposure time can be found at 

http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/pcarmf- http://cxc.harvard.edU/ciao/threads/streakextract/index.html#exposure 
11.7 http://cxc.harvard.edu/ciao/threads/spectra_hetgacis/ 
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The data analysis and model fitting throughout this pa- 
per were performed using XSPEC^^ version 12.6.0 (Arnaud 
1996) and, unless otherwise indicated, errors are everywhere 
reported at the la level of confidence. In this section (and 
throughout Appendices A and B), we fix the key input param- 
eters D, M and i at their fiducial values, which are given in 
Section 1. 

3.1. Seven Preliminary Models 

Our adopted model that is featured below, and upon which 
all of our results are based, was constructed by working in 
detail through a progression of seven preliminary models. We 
now briefly comment on these models, which are presented in 
full in Appendices A and B. 

Nonrelativistic models: The central component of our three 
nonrelativistic models, Models NR1-NR3, is the familiar 
accretion-disk model component DISKBB^^ (Mitsuda et al. 
1984; Makishima et al. 1986), which does not include any 
relativistic effects or the effects of spectral hardening, and 
which has an inappropriate boundary condition at the disk's 
inner edge (Zimmerman et al. 2005). We nevertheless em- 
ployed DISKBB as an exploratory tool because it has been 
widely used for decades, and it therefore allows us to com- 
pare our reduction/analysis results (for a spectrum of interest) 
to published results. Furthermore, this familiar model returns 
a direct and useful estimate of the temperature at the inner 
edge of the disk, which we use in order to make comparisons 
between Cygnus X-1 and other black hole binaries. The de- 
tails of this nonrelativistic analysis and the satisfying results 
obtained for Model NR3 - namely consistent values of the 
inner-disk radius and temperature for our three spectra - are 
presented in Appendix A. 

Relativistic models: Similarly, in addition to our adopted 
relativistic model (Section 3.2), in Appendix B we present our 
analysis and results for four relativistic models. Models Rl- 
R4, that are built around our fully relativistic accretion-disk 
model component KERRBB2, which we describe below. This 
component is a direct replacement for diskbb, returning two 
fit parameters, namely the spin and the mass accretion rate 
(instead of the temperature and radius of the inner disk). The 
four models progress sequentially in the sense that Model Rl 
is the most primitive and Model R4 is the most advanced. This 
sequence builds toward our adopted model. We have chosen 
to present our results for these preliminary relativistic models, 
in addition to those for our adopted model, because doing so 
demonstrates that our modeling of the critical thermal com- 
ponent, and the extreme spin it delivers for Cygnus X-1, are 
insensitive to the details of the analysis. 

3.2. Our Adopted Model 

The model we employ is a culmination of Models R1-R4 
in the sense that it is the most advanced and physically re- 
alistic model. The schematic sketch of the X-ray source in 
Figure 2 illustrates the various model components and their 
interplay. The structure of our adopted model, naming all the 
components that comprise it, is expressed as follows: 

CRABCOR * CONST * TEARS [SIMPLR (g) KERRBB2 
+KERRDISK + KERRCONV (g) (IREFLECT (g) SIMPLC)] 

XSPEC is available at http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 
For descriptions of the XSPEC models, see 
http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XspecModels.html 



As described in detail below, SIMPLR generates the power- 
law component using the seed photons supplied by the single 
thermal component KERRBB2, while the reflection compo- 
nent is likewise generated in turn by ireflect acting solely 
on the power-law component (i.e., ireflect does not act on 
the thermal component). Furthermore, the model fits for a sin- 
gle value of a* , which appears as the key fit parameter in three 
model components: KERRBB2, KERRDISK and KERRCONV. 

We now discuss in turn the model's three principal compo- 
nents - thermal, power-law and reflected - and their interrela- 
tionships. 

Thermal component: The centerpiece of our adopted model 
is our accretion-disk model KERRBB2, which includes all rel- 
ativistic effects, self-irradiation of the disk ("returning radia- 
tion") and limb darkening (Li et al. 2005). The effects of spec- 
tral hardening are incorporated into the basic model kerrbb 
via a pair of look-up tables for the hardening factor / cor- 
responding to two representative values of the viscosity pa- 
rameter: a = 0.01 and 0.1 (McClintock et al. 2006). Moti- 
vated by observational data obtained for dwarf novae (Smak 
1998, 1999) and soft X-ray transients (Dubus et al. 2001), 
and the results of global general relativistic magnetohydrody- 
namic (GRMHD) simulations, (Penna et al. 2010), through- 
out this work we adopt a = 0.1 as our fiducial value; mean- 
while, in Section 5.4 we examine the effects on our results 
of using a = 0.01 in place of a = 0.1. The entries in 
the look-up tables for / were computed using both kerrbb 
and a second relativistic disk model bhspec (Davis et al. 
2005; Davis & Hubeny 2006). We refer to the model ker- 
rbb plus this table, and the subroutine that reads it, as ker- 
RBB2 (McClintock et al. 2006). As noted above, the model 
KERRBB2 has just two fit parameters, namely the black hole 
spin parameter a* and the mass accretion rate M (or equiv- 
alently, a* and the Eddington- scaled bolometric luminosity, 

/ = I/boi(<^*, ^)/^Edd)- For the calculations reported in this 
paper, we included the effects of limb darkening and return- 
ing radiation. We set the torque at the inner boundary of the 
accretion disk to zero (as appropriate when D, M and i are 
held fixed), allowed the mass accretion rate to vary freely, and 
fitted directly for the spin parameter a*. 

Power-law component: The first term in the sum, 
SIMPLR(8) KERRBB 2, models the power-law component and 
the observed thermal component in combination. This domi- 
nant part of the spectrum (see Figure 3) is computed by con- 
volving KERRBB 2, which describes the seed photon distribu- 
tion (i.e., the thermal component prior to being scattered), 
with SIMPLR. The convolution model SIMPLR is a slightly 
modified version of SIMPL (Steiner et al. 2009b) that is appro- 
priate when including a separate and additive reflection com- 
ponent (Steiner et al. 2011). The parameters of SIMPLR (and 
SIMPL) are the power-law photon index F and the scattered 
fraction, /sc, which is the fraction of the seed photons that 
are scattered into the power-law tail. As used here, SIMPLR 
describes a corona that sends approximately half the scattered 
seed photons outward toward the observer and the remainder 
downward toward the disk, thereby generating the reflected 
component (see below). Thus, we assume that the power-law 
component illuminating the disk is the same as the component 
we observe. 

Reflected component: The second and third terms in the 
sum model the reprocessed emission from the disk that re- 
sults from its illumination by the power-law component. The 
model for the illuminating power-law component itself (the 
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term on the far right) is SIMPLC, which is equivalent to 
SlMPLR(g)KERRBB2 minus the unscattered thermal compo- 
nent (Steiner et al. 201 1). This power-law component is acted 
on by IREFLECT, which is a convolution reflection model with 
the same properties as its widely-used parent, the additive re- 
flection model PEXRiv (Magdziarz & Zdziarski 1995). Con- 
cerning IREFLECT, we (i) free the reflection scaling factor^^s 
while restricting its range to negative values, thereby com- 
puting only the reflected spectrum, while assuming that the 
corona is a thin slab that hugs the disk and emits isotropically; 
(ii) set the elemental-abundance switch to unity, which cor- 
responds to solar abundances, while allowing the iron abun- 
dance to be free; (iii) link the photon index to the value 
returned by SIMPLR; and (iv) fix the disk temperature at 
6.0 X 10^ K, the temperature Tin returned by diskbb (see Ap- 
pendix A). The model ireflect(8) SIMPLC returns a reflected 
spectrum that contains sharp absorption features and no emis- 
sion lines. To complete the model of the reflected component, 
we follow Brenneman & Reynolds (2006) and employ the line 
model KERRDISK and the convolution smearing model KER- 
RCONV, both of which treat a* as a free fit parameter^^. These 
models allow the emissivity indices to differ in the inner and 
outer regions of the disk. For simplicity, and because this 
parameter is unknown with values that vary widely from ap- 
plication to application, we use an unbroken emissivity pro- 
file with a single index q. We tie together all the common 
parameters of kerrdisk and kerrconv, including the two 
principal parameters, namely a* and q. 

The three multiplicative model components are, respec- 
tively, (1) CRAB COR, which corrects for detector effects (see 
Section 2); (2) CONST, which reconciles the calibration dif- 
ferences between the detectors (throughout the paper, we fix 
the normalization of the RXTEfPCUl detector and float the 
normalization of the ASCA GIS and Chandra HETG/ACIS 
detectors); and (3) TBABS^\ which models low-energy ab- 
sorption (Wilms et al. 2000). Concerning TBABS, throughout 
the paper we allow TVh to vary because the column density is 
well-determined by the data and TVh is known to vary by sev- 
eral percent for all three well- studied supergiant black-hole 
binaries, namely Cygnus X-1 (Hanke et al. 2009; also see Sec- 
tion 4), M33 X-7 (Liu et al. 2008) and LMC X-1 (Hanke et al. 
2010). 

4. RESULTS 

The fit results for our adopted model are given in Table 2, 
and the three fitted spectra together with their spectral com- 
ponents are shown in Figure 3. The fits are all good, with 
xl ranging from 1.17 to 1.28, and the results in Table 2 are 
in good agreement with those obtained using Model R1-R4 

The reflection scaling factor s in IREFLECT (and PEXRiv) is linked to the 
reflection constant parameter x in SIMPLR via the relation a::=l+|s| (where \s\ 
is the absolute value of s). In the limiting case of a: = 2, half of the scattered 
photons are redirected downward and illuminate the disc. In the limiting case 
of X = 1, none of the Compton- scattered photons strike the disc (Steiner 
et al. 2011). For the preliminary models described in the appendices that use 
the additive model PEXRiv, the scaling factor s is fixed at -1. 

We performed tests with a newer version of these models, relline and 
RELCONV (Dauser et al. 2010), which are interpolated on a finer grid, and 
found that our results presented in Section 4 (and elsewhere) are essentially 
unchanged. 

TBABS uses updated values for the photoelectric cross sections and ISM 
abundances and is an improved version of the familiar model PHABS. The 
choice of low-energy absorption model has a negligible affect on our results 
apart from the higher column density (~ 40%) that TBABS returns compared 
to PHABS. 



(see Appendix B and Tables 9-12). As in the case of Mod- 
els R1-R4, the spin parameter is high with a* > 0.99 for all 
three spectra. This is the principal result of this section. In 
the following section, we examine the robustness of this re- 
sult, and in Section 6 we fold in the uncertainties in the input 
parameters D, M and i and arrive at our final lower limit on 
a*. 

Another parameter of great interest is the scattering frac- 
tion, which measures the strength of the Compton compo- 
nent: fsc = 22.5 zb 0.6%, 30.5 zb 1.2% and 30.6 zb 0.6% for 
SPl, SP2 and SP3, respectively. These values are high com- 
pared to values characteristic of the thermal dominant state 
(/sc ~ 5%; Steiner et al. 2009b), which is the most favorable 
state for the measurement of black hole spin. These large val- 
ues of /sc? as well as the uniform value of the power-law in- 
dex (F ~ 2.5), imply that for all three observations the source 
was in the steep power-law state (Remillard & McClintock 
2006). Steiner et al. (2009b) have shown that one can obtain 
reliable values of the inner disk radius (and hence spin) in the 
SPL state for fsc ^ 25% (see their Figure 1). 

The luminosity of the disk component, L/L^dd ~ 0.02, 
is quite low and easily meets our data selection criterion, 
L/LEdd < 0.3 (McClintock et al. 2006). Correspondingly, 
the disk is geometrically very thin at all radii (the aspect ra- 
tio h/r < 0.05; see Penna et al. 2010; Kulkarni et al. 2011). 
At the same time, the luminosity is sufficiently high that the 
spectral hardening factor / is well-determined by the data 
(/ ^ 1.6). 

Interestingly, for the pair of Chandra observations that are 
separated by just two days, the values of the fit parameters 
are quite similar (Table 2) with the notable exception of the 
column density: A^h = 0.768 ± 0.024 cm-^ for SP2 and 
A^H = 0.687 zb 0.010 cm-2 for SP3. We tested whether one 
can achieve a good fit with a single value of A'h, by fitting 
the Chandra/RXTE spectra jointly, linking only the parameter 
A^H- This linking increased the total chi-square by 41, which 
corresponds to an F-test probability of 3.0 x 10~^. This low 
probability implies a significant change in A'h, which is likely 
due to absorption in the wind of the supergiant companion. 
Indeed, one expects such absorption in a wind to be larger 
for SP2 at orbital phase 0.24 than for SP3 at phase 0.61, as 
observed (Table 1; Balucihska-Church et al. 2000). 

5. ROBUSTNESS OF SPIN ESTIMATES 

In this section we consider a number of factors that might 
affect our results and find that none of them is significant. 
We consider the effects of: (1) excluding the Fe Ka line and 
edge features from the fits; (2) including HEXTE data and 
extending the fits to 150 keV; (3) using REFLIONX to model 
the reflection component; and (4) substituting a = 0.01 for 
our fiducial value of a = 0.1 and varying the metallicity. As 
shown below, factors (1) and (2) have negligible effects on 
our results, and (3) and (4) have slight upward effects on a*, 
which implies that the extreme values reported in Table 2 are 
conservative lower limits. Finally (Section 5.5), we explore 
the effects on our results of artificially relaxing the spin pa- 
rameter away from its extreme, limiting value by varying the 
parameters of KERRBB2; we find that the input parameters D, 
M, and i have to be driven far from their fiducial values in 
order to obtain a spin value as low as a* =0.9. 

5.1. Effect of iron line and edges 

For all three spectra, we refitted the data excluding 
the energy range 5.0-10.0 keV while omitting the compo- 
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nent KERRDISK. This excised energy range contains the 
relativistically-broadened Fe Ka Hne and edge, as well as a 
significant residual feature near 9 keV (Figure 3). For all three 
spectra, we find that our results are essentially unchanged, 
apart from small shifts in the parameters of the reflection com- 
ponent. The results for spectrum SPl (only) are shown in 
Table 3 (Case 2) where they are compared to our standard 
results (Case 1). Thus, we find that the values of the spin 
parameter returned by the fits are completely determined by 
the temperature and luminosity of the thermal component and 
are unaffected by the presence of the line. This conclusion 
is reasonable given that the line is a minor feature with an 
equivalent width of 0.15 keV, or only ^ 0.015 keV relative 
to the continuum at the peak of the thermal component (see 
Figure 3). 

5.2. Effect of extending the bandwidth to 150 keV 

In Section 2, we asserted that the coverage of the PCA to 
45 keV was sufficient to adequately constrain the power-law 
and reflection components. We now demonstrate that this is 
true by refitting the data for SPl (only) while including the 
HEXTE data spanning the energy range 20 keV to 150 keV 
(where the source counts are negligible). For the HEXTE 
data, we do not a add a systematic error to the count rates. 
However, we do correct the detector response to the Crab 
spectrum of Toor & Seward (see Section 2) using the value 
of the Crab's photon index as measured using HEXTE^^: 
F = 1.93 zb 0.003. The slope correction is AFts = -0.17. 
The results obtained using the HEXTE data, which are given 
in Table 3 (Case 3), are almost identical to our standard results 
for SPl (Case 1). This is not surprising because the PCA cov- 
erage to 45 keV is more than adequate to determine the slope 
of the power-law component, and the reflection component is 
quite weak and dying rapidly at 45 keV (Figure 3). 

5.3. Effect of using a different reflection model 

We replace IREFLECT(g)SIMPLC+KERRDISK with RE- 
FLIONX^-^ (Ross & Fabian 2005), which is widely used in 
measuring spin via the Fe Ka line. The merit of this alter- 
native reflection model is that it self-consistently calculates 
the line feature and the reflection component, whereas IRE- 
FLECT models only the absorption edges. The downside of 
REFLIONX is its description of the seed photon distribution as 
a simple power-law, which unphysically diverges at low en- 
ergies (for further comparison of the two models, see Steiner 
et al. 2011). We include relativistic blurring by convolving 
REFLIONX with KERRCONV. Compared to the results given 
in Table 2, the values of the parameters returned by the fits 
are generally different, although reasonable, and the fits are 
somewhat poorer (x^ =1.35-1.45). Importantly, the effects 
on the spin parameter are very small: The value of a* in- 
creases slightly for SPl (0.9985 to 0.9999) and is unaffected 
for SP2 and SP3, which are at their maximum values. 

5.4. Effect of Varying the Viscosity Parameter and Metallicity 

We refitted the spectra using a = 0.01 in place of our fidu- 
cial value, a = 0.1. Again, the spin parameter for SPl in- 
creases sHghtly (0.9985 to 0.9988), while the values for SP2 

22 http://web.mit.edu/iachec/IACHEC_2_talks/IACHEC_II_suchy.pdf 

23 The model can be downloaded at 
http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/reflion.html 



and SP3 remain unchanged. Concerning possible metallic- 
ity effects, we do not have nonsolar-metallicity table mod- 
els for computing spectral hardening for such extreme val- 
ues of spin. However, for three sources (M33 X-7, Liu 
et al. 2008; LMC X-1, Gou et al. 2009 and A0620-00, Gou 
et al. 2010), we have found that substantial changes in metal- 
licity produce very small changes in a*. For example, re- 
ducing the metallicity from solar to a tenth solar decreases 
the spin parameter of a high-spin black hole like LMC X-1 
(a* = 0.938 zb 0.020) by only Aa* = 0.001. The effect of 
this same change in metallicity for the slowly- spinning black 
hole A0620-00 (a* = 0.135 =b 0.029) is larger, but it is still 
quite small: Aa* = 0.014 . (The small errors on a* quoted 
here come directly from fitting the X-ray spectra, and they 
do not take into account the uncertainties in D, M and i, 
which dominate the error budget). We note that the super- 
solar iron abundances implied by our fits using the reflection 
models, ireflect (Table 2) and pexriv (Appendix B), sug- 
gest that the metallicity of the Cygnus X-1 system is possibly 
enhanced. Accounting for the supersolar abundances will re- 
sult in a slightly increased estimate of a*, so our conclusions 
are robust to enhanced metallicity. 

5.5. Relaxing the Spin Parameter 

We now describe three technical exercises that artificially 
examine the effects of varying three parameters of KERRBB2, 
namely, a^, D and the normalization constant A/'k (which we 
have elsewhere fixed to unity, as is appropriate when D, M 
and i are specified). Our motivation is to examine the conse- 
quences of relaxing the spin parameter away from the extreme 
values returned by the fits (Table 2). 

First, for spectrum SPl only, we leave A'k fixed at unity and 
allow the distance to vary (keeping M and i fixed at their fidu- 
cial values). For our parallax distance ofD = 1.86 kpc, we of 
course have our standard result, a* = 0.9985 (Table 2). We 
now arbitrarily and successively fix the spin parameter at two 
lower values, a* = 0.95 and a* = 0.90, and refit spectrum 
SPl. The best-fit values of D are then substantially greater 
than our measured value of 1.86 kpc: 2.43 ± 0.03 kpc and 
2.78 =b 0.02 kpc, respectively. Meanwhile, the corresponding 
values of reduced chi-square are respectively 1.34 and 1.40, 
which are significantly greater than our standard value of 1.24 
(Table 2). Thus, the best fit is achieved for the observed and 
extreme value of spin. 

Secondly and alternatively, we obtain similar results by 
varying N^^ while leaving the distance fixed at its fiducial 
value of 1.86 kpc. For the same pair of forced values of 
the spin parameter given above (0.95 and 0.90), we find for 
spectrum SPl that the fitted values of A/'k are about half the 
standard value of unity: 0.65 ± 0.01 and 0.45 ± 0.01, re- 
spectively (where a smaller value corresponds to a weaker 
thermal component). Meanwhile, we find values of reduced 
chi-square that are very similar to those given in the previous 
example, with values that increase as the spin parameter de- 
creases. That is, we again find that the best fit is achieved for 
the observed and extreme value of spin. 

In a final experiment, we assess the consistency of our spin 
values for the three spectra, which can not be adequately 
judged by considering the extreme values given in Table 2. 
We make this assessment by artificially setting A/'k = 0.5 in 
our adopted model and refitting the three spectra. In this way, 
we find best-fit spin values that are in reasonable agreement: 
a* = 0.937t^:^^} for SPl, a* = 0.934l^;^^J for SP2, and 



The Extreme Spin of Cygnus X- 1 



7 



a* = 0.9531^;^^^ for SP3. 

6. COMPREHENSIVE ERROR ANALYSIS 

In all previous work on measuring spin, we have found 
that the statistical uncertainty in the estimated spin parame- 
ter due to the X-ray data analysis is small compared to that 
due to observational uncertainties in D, M and i. Based on 
our GRMHD simulations of thin disks (Shafee et al. 2008; 
Penna et al. 2010), we have likewise found that the uncer- 
tainties in these three input parameters dominate over the er- 
rors resulting from our use of the analytic Novikov-Thorne 
model, which assumes a razor-thin disk. In the case of Cygnus 
X-1, these model errors are especially small because of the 
extreme spin of the black hole and the low luminosity of 
the disk. Spin estimates obtained by fitting mock spectra 
of simulated GRMHD disks indicate that for an inclination 
i = 30 deg, which is very near the inclination of Cygnus X-1 
(i = 27.1 deg), the Novikov-Thorne thin-disk model overesti- 
mates the spin parameter by only Aa* = 0.007 and 0.005 for 
spins of 0.90 and 0.98, respectively (see Table 1 in Kulkarni 
et al. 2011). Furthermore, these errors are significantly over- 
estimated because they were computed for disks that are far 
more luminous (L/LEdd ^ 0.5), and hence thicker, than that 
of Cygnus X-1 (L/i^Edd ^ 0.02)^"^. In our error analysis, we 
neglect this small model error. 

With the exception of our recent spin measurement of XTE 
J1550-564 (Steiner et al. 2011), in earlier work we have ne- 
glected the uncertainty in the luminosity due to the 10% 
uncertainty in the flux of the Crab (Toor & Seward 1974), an 
error that uniformly shifts all of our spin measurements either 
up or down. For the spin of Cygnus X-1, the effect of the 
uncertainty in the absolute flux calibration is very compara- 
ble to the 6% uncertainty in D (which is equivalent to a 12% 
uncertainty in the measurement of flux). We have therefore in- 
cluded in our error budget the 10% uncertainty in flux (which 
we approximate as an uncertainty in the distance of 0.1 kpc) 
by simply combining the distance and flux-calibration errors 
in quadrature, thereby inflating the actual 0.120 kpc distance 
uncertainty to 0.156 kpc. Thus, the final error we report for 
a* includes the uncertainties in D, M, i and the uncertainty 
in the absolute flux calibration. Taken together, these four 
sources of uncertainty totally dominate the error budget. 

In order to determine the error in a* due to the combined 
uncertainties in D, M and i, we performed Monte Carlo 
simulations using the Odyssey computing cluster at Harvard 
University. The latter two parameters are not independent. 
They are related through the expression for the mass func- 
tion: /(M) = M^^ixt'i/^M + Mopt)^ = 0.263 ± 0.004 M©, 
where Mopt = 19.16 ± 1.90 Mq is the mass of the secondary 
star and the value of the mass function was evaluated using a 
K- velocity of 76.79 =b 0.41 km and an orbital period of 
P = 5.599829 days (Orosz et al. 2011). 

In the analysis, we assumed that the value of the mass func- 
tion, the inclination, and the mass of the secondary are nor- 
mally and independently distributed, and we computed the 
mass of the black hole using the values of these quantities, 
which are given above. We conservatively fixed the viscos- 
ity parameter at its baseline value, a = 0.1 (using a = 0.01 
increases a*; see Section 5.4). Specifically, we (i) generated 
9000 parameter sets for D, i, Mopt, and f{M) \ (ii) solved for 
M for a given triplet of values of i, Mopt and /(M) ; (iii) com- 
puted for each set the look-up table for the spectral hardening 

It is computationally very challenging to simulate thinner disks. 



factor / using the model BHSPEC; and (iv) obtained a* by 
fitting our adopted model to the spectra. The final histogram 
distributions for our three spectra are shown in Figure 4. Con- 
sulting these histograms, we see that the spin estimate is low- 
est for SPl, and we conservatively base our final result on this 
spectrum, thereby concluding that a* > 0.95 a? the 3cr level 
of confidence (see Figure 4 and the bottom line of Table 2). 

We note the following two caveats: First, the power- 
law component is strong relative to the thermal compo- 
nent, which decreases the reliability of the continuum-fitting 
method (Steiner et al. 2010, 2011). Second, as we discuss 
in Section 7.4, the continuum- fitting method assumes that the 
spin vector of the black hole is aligned with the orbital angular 
momentum vector. 

7. DISCUSSION 

Our wide-ranging discussion covers seven topics: (1) We 
first confront the challenge posed by a grossly discrepant spin 
result that was obtained using the Fe-line method (while not- 
ing a concordant Fe-line result that appeared very recently), 
and (2) a second discrepant result obtained using a QPO 
method. (3) We next consider a disfavored dynamical model 
that implies a less extreme value of spin. (4) The fundamen- 
tal assumption of the continuum-fitting method, namely the 
alignment of the spin and orbital vectors, is then considered. 
(5) The extreme spin of the black hole is reconciled with the 
low temperature of its accretion disk. (6) We find no evidence 
in our Chandra HETG spectra for warm absorbing gas, and 
we show that, if it were present, its effects on our spin esti- 
mates would be negligible. (7) Finally, we discuss some con- 
sequences of the extreme spin of Cygnus X-1, and we argue 
that this spin is chiefly natal in origin. 

7.1. Measurement of Spin via the Fe Ka Method 

Based on an analysis of the Fe line profile. Miller et al. 
(2005) found marginal evidence for high spin using the same 
ASCA spectrum of Cygnus X-1 that we use. Subsequently, 
however. Miller et al. (2009) reported a spin of a* = 0.05 i 
0.01 based on an analysis that combines Fe Ko^/reflection 
models and continuum models, including kerrbb, a rela- 
tivistic disk model similar to the one we use. That is, the au- 
thors jointly applied the continuum- fitting and Fe Ka methods 
and fitted simultaneously for a single value of a* . The data an- 
alyzed in this case were a pair of 0.5 ks XMM-Newton/fE?\C- 
pn spectra obtained in the "burst" timing mode. 

While it is not possible for us to account in detail for the 
gross difference between the near-zero spin reported by Miller 
et al. (2009) and our near-extreme value, we note the fol- 
lowing: First, the round values of D, M, and i that Miller 
et al. used as input to kerrbb differ significantly from ours, 
and these differences all serve to drive the spin down; e.g., 
using their values of these parameters and our spectrum SPl, 
we find a* = 0.74 ±0.01. Secondly, and in addition. Miller 
et al. fitted for the kerrbb normalization constant, obtain- 
ing A^K = 0.31, whereas the standard procedure, which we 
follow, is to set A/'k to unity when D, M, and i are fixed. Fit- 
ting SPl using TVk = 0.31 (and the values of D, M, and i 
adopted by Miller et al.), we find that the spin drops from 
a* = 0.78 to the retrograde value a* = —0.53 (and the fit is 
poor, =2.01.) 

The near-zero spin result of Miller et al. (2009) is further 
called into question because no data above 10 keV were used. 
This lack of high energy coverage, in the presence of a strong 
Compton component, seriously compromises results obtained 
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using the continuum- fitting method, as we stress at the outset 
of Section 2 (and as can be deduced by an examination of Fig- 
ure 3). Likewise, results obtained using the Fe Ko^/reflection 
method are compromised by the failure to observe the Comp- 
ton reflection hump around 20-30 keV (e.g., see Larsson et al. 
2008). 

Very recently, after posting our paper on the astro-ph 
archive and while revising it for resubmission, a paper ap- 
peared reporting another estimate of the spin of Cygnus X-1 
via an analysis of the Fe Ka profile (Duro et al. 2011). As- 
suming that the emissivity profile of the disk can be described 
by a single power-law with index 3.0, these authors conclude 
that "the black hole is close to rotating maximally," which is 
in agreement with our result. We note that Duro et al. do not 
discuss Miller et al. (2005, 2009) or mention the near-zero 
spin result reported in the latter paper. 

7.2. Measurement of Spin via a QPO Model 

Based on an analysis of low-frequency (0.01-25 Hz) quasi- 
periodic oscillations (QPOs), Axelsson et al. (2005) obtained 
a spin for Cygnus X-1 of a* = 0.49 =b 0.01 for M = SM©. 
Their result is based on the relativistic precession model of 
Stella et al. (1999). In this model, the QPO is produced as a 
result of emission from an orbiting bright spot that is undergo- 
ing relativistic nodal and periastron precessions in a slightly 
tilted and eccentric orbit. The spin parameter is predicted to 
vary with mass as a* oc M~^/^ (see Eqn. 4 in Axelsson et al. 
2005), and therefore the corrected value of spin is a* = 0.43 
for our adopted black hole mass M = 14.8 Mq. The large 
discrepancy between this moderate value of spin and the ex- 
treme value we find may be a consequence of a fundamental 
assumption of their model, namely, that the black hole is ro- 
tating slowly (a* <^ 1). On the other hand, the precession 
model, with its assumption of geodesic motion, may not ap- 
ply in this instance. 

7.3. An Alternative Dynamical Model 

In Orosz et al. (201 1), results are presented for four dynam- 
ical models. Models A — V. We disregard Models A and 
S, which assume a circular orbit, because these models give 
poor fits compared to the eccentric orbit models (Ax^ > 50) 
and because there is clear evidence that the orbit is eccentric. 
Throughout this paper, we have used M = 14.8±1.0 Mq and 
i = 27.1 ± 0.8 deg from Model V, an asynchronous model 
with a rotational frequency for the 0-star that is 40% greater 
than the orbital frequency. As an alternative to Model V, we 
now consider Model C, which assumes synchronous rotation. 
Model C gives a poorer fit to the data (Ax^ ~ 13), and it re- 
sults in disharmony between the light-curve and velocity data 
on the one hand, and the radius and rotational velocity of the 
0-star on the other (see Table 1 in Orosz et al. 201 1). Because 
of this disharmony, the uncertainties in M and i for Model C 
are significantly larger than those for our favored model, al- 
though the central values of these parameters differ only mod- 
estly: M = 15.8 ± 1.8 Mq and i = 28.5 ± 2.2 deg. Using 
these values for Model C and spectrum SPl (which gives the 
lowest value of spin), we repeated our error analysis (Section 
6) and obtained the following 3 a lower limit on the spin pa- 
rameter: a* > 0.92. 

7.4. Alignment of Spin and Orbital Angular Momentum 

There remains one uncertainty that calls into question the 
reliability of the continuum-fitting method, namely, whether 



the inner X-ray emitting portion of the disk (which will align 
with the black hole's spin axis) is aligned with the binary or- 
bital plane. For a discussion of this question, see Section 2.2 
in Li et al. (2009). As an extension of this discussion, we 
note the following: First, recent population synthesis stud- 
ies predict that the majority of systems will have rather small 
(< 10 deg) misalignment angles (Fragos et al. 2010). Sec- 
ond, in the case of Cygnus X-1, there is reason to believe 
that the misalignment angle is especially small because of 
the binary system's low peculiar velocity, which indicates that 
the system did not experience a large "kick" when the black 
hole formed (Mirabel & Rodrigues 2003; Reid et al. 2011). 
As demonstrated in Figure 5, even if there exists a misalign- 
ment angle as large as (e.g.) 16 degrees, the spin value is still 
>0.95. 

7.5. Low Disk Temperature Compatible with Extreme Spin 

One might expect that the accretion disk of a fast- spinning 
black hole like Cygnus X-1 would be hot (Tin >1 keV) be- 
cause its inner edge is relatively close to the event horizon 
and deep in the gravitational potential well. Two principal 
factors are responsible for the depressed disk temperature 
(Tin ^0.5 keV; Table 7). The first of these is the low rate 
of mass accretion through the disk, which is manifested by its 
low luminosity, only ^ 2% of Eddington (Table 2). The sec- 
ond is that relativistic effects (e.g., beaming and light bend- 
ing) are muted for low disk inclinations: Using our fiducial 
values of D, M, and i, and fixed values of the fit parameters 
taken from Table 2, we used KERRBB2 to simulate a pair of 
spectra, one for the inclination of Cygnus X-1, i = 27.1 deg, 
and the other for i = 66.0 deg (which is the inclination of 
GRS 1915+105; McClintock et al. 2006). The luminosities 
of the two spectra are the same, L/LEdd ^ 0.018, while the 
disk temperatures are Tin = 0.45 keV and 0.86 keV, respec- 
tively. Thus, in this comparison, the temperature of the low- 
inclination disk is depressed by nearly a factor of two because 
the relativistic effects are weak. 

7.6. Effects of a Warm Absorber 

Although careful modeling of warm absorbers (WAs) is 
usually crucial in determining the spins of supermassive black 
holes via the Fe Ka method (e.g., Brenneman & Reynolds 
2006), we find that the effects of WAs are unimportant in 
estimating the spin of Cygnus X-1 via the continuum-fitting 
method. We examined all of the available soft- state Chan- 
dra HETG spectra of Cygnus X-1 (ObsIDs 2741, 2742, 2743, 
12313/SP2, and 12314/SP3) Sit E < 1 keV, and we find no 
evidence for the blend of absorption lines near 0.76 keV (due, 
e.g., to Fe X-Fe XV, O VII and O VIII), which is a telltale 
signature of a WA. (The lines WAs produce above 1 keV are 
mostly discrete and weaker and therefore have an accordingly 
much smaller effect on the continuum shape and our spin es- 
timates.) 

Our preliminary analysis of the low-energy portion (E < 
0.8 keV) of a 50 ks hard-state Chandra HETG spectrum 
(Hanke et al. 2009) did reveal the presence of a single WA. 
Its two most relevant parameters are its column density, A^h = 
7.2 X 10^° cm~^, and ionization parameter, ^ = 44. (Its turbu- 
lent broadening and redshift are respectively v = 47.6 km s~^ 
and z = —0.00086; we assume solar abundances.) This ten- 
uous WA does not affect the continuum shape, but it does 
produce line blends, the strongest of which is centered at 
0.76 keV. Although such a WA is not present in soft-state 
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spectra, we nevertheless tested its effects by reanalyzing SPl, 
SP2 and SP3 by introducing the additional model compo- 
nent WARMABS, fixing its parameters to the values given 
above; we obtained values of spin that are essentially identi- 
cal to those given in Table 2. In one further test, we likewise 
reanalyzed our three spectra including a second thicker and 
hotter WA that generates a complex of absorption lines be- 
tween 0.9 and 1.0 keV (A^h = 1.0 x 10^^ cm-^; ^ = 300; 
v = 100 km s-^ z = -0.00086). Again, the effects on our 
spin estimates are negligible. 

7.7. Consequences of Extreme Spin and its Origin 

The spin of Cygnus X-1 is extraordinarily high, placing it 
in the company of the microquasar GRS 1915+105, estimated 
using the continuum-fitting method (McClintock et al. 2006; 
Blum et al. 2009), and the supermassive black hole in MCG- 
6-30-15, estimated using the Fe-line method (Brenneman & 
Reynolds 2006). The spins of both of these black holes are 
reported to exceed a* = 0.98. As Penrose (1969) first demon- 
strated, the enormous "flywheel" energy of a fast- spinning 
black hole can in principle be tapped. For Cygnus X-1, the 
potentially tappable energy is >2.8 Mqc^ = 5.0 x 10^^ ergs 
(Christodoulou & Ruffini 1971); in comparison, the energy 
radiated by the Sun over its entire ~ 10 billion-year lifetime 
is < 0.001 MqC^. It has been widely suggested that spin- 
down energy powers the relativistic jets observed for at least 
some quasars and microquasars (Blandford & Znajek 1977), 
such as GRS 1915+105 (Mirabel & Rodriguez 1994). 

As the spin a* of a black hole approaches unity, the radius 
of its ISCO approaches the radius of the event horizon R^u 
(Bardeen et al. 1972). For Cygnus X-1, with spin a* > 0.95 
and M = 14.8 M©, i^isco < 42 km while Reu < 29 km. 
The Keplerian velocity of the gas at the ISCO is approxi- 
mately half the speed of light, and its orbital frequency is 

> 598 Hz. Meanwhile, the frequency of rotation of the black 
hole itself, which is the spin frequency of space-time at its 
horizon, is > 790 Hz (the maximal spin frequency for a* = 1 
is 1091 Hz). 

What is the origin of the spin of Cygnus X- 1 ? Was the black 
hole born with its present spin, or was it torqued up gradually 
by the gas it has accreted over its lifetime? To achieve a spin 
of a* > 0.95 via disk accretion, an initially nonspinning black 
hole must accrete > 7.3 Mq from its donor (Bardeen 1970; 
King & Kolb 1999) in becoming the M = 14.8 M© that 
we observe today. However, to transfer this much mass even 
at the maximum (Eddington-limited) accretion rate requires 

> 31 million years^^, whereas the age of the binary system is 
between 4.8 and 7.6 million years^^ (Wong et al. 201 1). (Even 
for a* > 0.92 obtained for our less probable model, the accre- 
tion timescale is > 25 million years.) Thus, it appears that the 
spin of Cygnus X-1 must be chiefly natal (also, see Axelsson 
et al. 201 1 ; Wong et al. 201 1), although possibly the high spin 
could be achieved during a short-lived evolutionary phase of 
hypercritical mass accretion (Moreno Mendez 201 1). 

The corresponding estimates of spin-up time we reported for both 
M33 X-1 (Liu et al. 2008, 2010) and LMC X-1 (Gou et al. 2009) are in- 
correct because for the efficiency r] (of converting mass to radiant energy) we 
used a constant value, r] = 1. However, gradually increases as the black 
hole spins up, starting at a* =0 with 77 = 0.06 to r] = 0.13 (a* = 0.84) 
for M33 X-7 and r] = 0.17 (a* = 0.92) for LMC X-1. For M33 X-7 and 
LMC X-1, the correct limits on the spin-up timescales are > 17 and > 25 
million years, respectively, while the respective ages of the systems are < 3 
and < 5 million years. 

A 7.6 million-year-old system accreting at the maximum rate would 
only achieve a spin of 0.59. 



8. CONCLUSION 

Based on our favored dynamical model, we find an ex- 
treme value of spin for the black hole primary in Cygnus X-1 : 
a* > 0.95 at the 3a level of confidence. For a less prob- 
able (synchronous) dynamical model, the spin is still high: 
a* > 0.92 (3a). For both of these strong limits, we include 
the customary uncertainties in the input parameters D (6%), 
M (7%) and i (±0.8 deg), and we also include the uncertainty 
in the calibration of the absolute flux via the Crab (10%). 
These four sources of uncertainty totally dominate the error 
budget. 

Our measurement of spin is determined solely by the prop- 
erties of the thermal component and is unaffected by the pres- 
ence of the relatively faint Fe Ka line. Nevertheless, we have 
modeled this relativistically-broadened line feature carefully 
in order to achieve good fits over the full range of energies 
we consider, which in our analyses variously extends from 
0.5 keV to 150 keV. The extreme spin we find for this black 
hole is based on an analysis of three broadband spectra that 
are each capable of constraining the soft thermal component, 
the hard Compton component, and the reflected component. 
By considering several different models and performing a 
number of tests on our results, we have demonstrated that the 
extreme spin we find is insensitive to the details of our analy- 
sis. 
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1996 1998 2000 2002 2004 2006 2008 2010 

Years 

Figure 1. Fourteen-year record for Cygnus X-1 of X-ray intensity relative to the Crab nebula (top) and spectral hardness SH (bottom). The hardness is defined 
as the ratio of counts detected in a hard X-ray band (5-12 keV) to those detected in a soft (1.5-5 keV) band. We consider data suitable for the measurement of 
spin only when the spectral hardness is below the dashed line (SH < 0.7), which is an empirical choice. Shown plotted as red stars are the intensity and hardness 
of the source for each of the three selected observations, SPl in 1996, and SP2 and SP3 in 2010. While a useful diagnostic for the purposes of data selection, the 
RXTE/ ASM sky survey data shown here are unsuitable for the measurement of spin. 
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Power-Law Component 




Figure 2. Schematic sketch of the X-ray source (adapted from a sketch provided by R. Reis). The accretion disk (pink) is truncated at the ISCO, leaving a dark 
gap between the disk's inner edge and the black hole's event horizon (black). Shown hovering above the optically-thick disk is its tenuous scattering corona 
(yellow). As indicated by the arrows, the disk supplies the thermal component of emission, which is Compton scattered into a power-law component by hot 
electrons in the corona. Approximately half of this latter component illuminates the disk, thereby generating the reflected component. 




Figure 3. (Top) The upper envelope in each of these spectra shows the data (RXTE in blue, and ASCA or Chandra in black) and the best-fit total relativistic 
model for the case of our adopted model. Each total model spectrum is shown decomposed into thermal and power-law components, and a reflection component, 
which is comprised of a continuum component plus the F Ka line feature. (The color assignments correspond to those used in Figure 2.) The low-energy X-ray 
absorption component is evident at energies <1 keV. Note in all three spectra the dominance at low energies of the key thermal component. (Bottom) Ratio of 
the data to the model showing deviations between the two. 




Figure 4. Histograms of a* for 9000 parameter sets (per spectrum) resulting from the Monte Carlo analysis. The lower limits given are at the 3cr (99.7%) level 
of confidence. 
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Figure 5. The effect on the spin parameter of varying the incHnation angle i (for SPl only) for fixed values of our adopted parameters: D = 1.86 kpc and 
M = 14.8 Mq. The best-fit value of inclination (i = 27.1 deg) is indicated by the dashed line. While the spin parameter in this figure varies modestly from 
0.93 to unity, the ISCO radius - the quantity that we actually measure - can be seen to change by a factor of 2 (Bardeen et al. 1972). As an extreme example, 
as a* increases from 0.98 to 0.99, the fractional change in the ISCO radius is 10 times the fractional change in the spin parameter. This "saturation" of the spin 
parameter near unity is the reason that the statistical uncertainties given in Table 2 are so small. 
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Table 1 

Journal of the observations that yielded spectra SPl, SP2 and SP3' 



No. 


Mission 


Detector 


E1-E2 (keV) 


UT 


^exp (sec) 


I(Crab) 


SH 


4> 


SPl 
SP2 
SP3 


ASCA & RXTE 
Chandra & RXTE 
Chandra & RXTE 


GIS & PCA 
ACIS(TE)/HETG & PCA 
ACIS(CC)/HETG & PCA 


0.7-8.0 &2.5-45.0 
0.5-10.0 & 2.5-45.0 
0.5-10.0 & 2.5-45.0 


1996-05-30 06:43:16 & 07:51:29 
2010-07-22 16:21:22 & 17:04:01 
2010-07-24 17:21:43 & 17:43:00 


2547 & 2240 
2146 & 3808 
900 & 3904 


0.80 
1.31 
1.01 


0.64 
0.54 
0.61 


0.74 
0.24 
0.61 



^ For each pair of missions listed, in columns 3-6 we give respectively the following information: names of the detectors employed, bandwidths used in the 
analysis, UT start times of the observations and effective exposure times. The source intensity and spectral hardness (SH), which are defined and plotted in 
Figure 1, are given in the two columns that follow. The final column gives the orbital phase of the binary system defined as the time of supergiant superior 
conjunction (black hole beyond O-star), which occurred on hehocentric Juhan Day 2441874.71 (Brocksopp et al. 1999). 



Table 2 

Fit results for our adopted model' 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 


KERRBB2 


a* 


U.»yOO_Q QQQg 


n QQQQ + O b 

U.»»»»_Q QQ29 


QQQQ + b 


2 


KERRBB2 


M 


0.115 ±0.004 


0.139 ±0.029 


0.122 ±0.011 


3 


const 




1.000 ±0.002 


1.031 ±0.013 


0.971 ±0.004 


4 


TBABS 




0.705 ± 0.006 


0.768 ± 0.024 


0.687 ±0.010 


5 


SIMPLR 


r 


2.282 ±0.010 


2.499 ±0.012 


2.549 ±0.011 


6 


SIMPLR 


/sc 


0.225 ±0.006 


0.305 ±0.012 


0.306 ± 0.006 


7 


KERRDISK 




6.56 ±0.09 


6.44 ± 0.05 


6.49 ±0.04 


8 


KERRDISK 


q. 


2.82 ±0.02 


2.28 ±0.06 


2.29 ±0.07 


9 


KERRDISK 




0.015 ±0.001 


0.032 ± 0.002 


0.024 ± 0.001 


10 


KERRDISK 


EW 


0.154 


0.158 


0.159 


11 


IREFLECT 




5.34 ±0.15 


2.97 ±0.12 


3.63 ±0.14 


12 


IREFLECT 


S 


0.98 ± 0.04 


0.80 ± 0.03 


0.93 ± 0.03 


13 


IREFLECT 




153.1 ±15.7 


71.2 ± 10.0 


62.0 ±8.4 


14 
15 




xi 

f 


1.24(561/454) 
1.610 


1.28(371/290) 
1.612 


1.18(723/614) 
1.621 


16 






0.018 


0.026 


0.023 


17 


Adopted ^ 


a* 


9985+^-^^^^ 


0.9999^^.0050 


0.9999^^.0116 


^ For the model components given, 


the parameters from top to bottom are: 


(1) spin parameter; (2) mass accretion rate in units of 10^^ 


g s~^; (3) detector 



normalization constant relative to RXTE PCU2; (4) hydrogen column density in units of cm~^; (5) photon power-law index F; (6) scattering fraction /sc; 
(7) central line energy in keV; (8) emissivity index q\ (9) line flux in units of photons cm~^ s~^; (10) equivalent width of line in keV; (11) iron abundance 
relative to solar; (12) reflection scaling factor s; (13) ionization parameter (14) Reduced chi-square, total chi-square and degrees of freedom, respectively; (15) 
spectral hardening factor /; and (16) Eddington-scaled disk luminosity, where Lpdd ~ 1-9 x 10^^ erg for Cygnus X-1. Unless otherwise indicated, the 
uncertainties quoted here and throughout the paper are at the Icr level of confidence. 

^ The physical maximum value of the spin parameter a* is unity and for the XSPEC model kerrbb2 it is 0.9999. The very small errors quoted here are purely 
the uncertainties due to counting statistics, which result from fitting our adopted model to the X-ray data. 

^ Final adopted values for the spin parameter and their uncertainties. The Icr uncertainties are calculated based on the 3cr lower limits on a* shown in Figure 4. 
These results fold in the uncertainties in D, M,i and the absolute flux calibration via our Monte Carlo analysis (see Section 6). 



Table 3 

Effects of excluding the Fe line and increasing bandwidth (SPl only)^ 



Number 


Model 


Parameter 


Case 1 


Case 2 


Case 3 


1 


KERRBB2 


a* 


9985+°-°°°^ 




u.yyo4_Q QQQ3 


2 


KERRBB2 


M 


0.115 ±0.004 


0.115 ±0.007 


0.115 ±0.003 


3 


const 




1.000 ±0.002 


0.999 ± 0.002 


1.000 ±0.002 


4 


TBABS 


Nn 


0.705 ± 0.006 


0.708 ± 0.009 


0.706 ± 0.005 


5 


SIMPLR 


F 


2.282 ±0.010 


2.288 ±0.013 


2.287 ±0.008 


6 


SIMPLR 


/sc 


0.225 ± 0.006 


0.221 ±0.010 


0.221 ± 0.005 


7 


KERRDISK 


^L 


6.56 ±0.09 


0.00 ±0.00 


6.56 ±0.05 


8 


KERRDISK 


q 


2.82 ±0.02 


2.82 ±0.23 


2.81 ±0.02 


9 


KERRDISK 


A^L 


0.015 ±0.001 


0.000 ± 0.000 


0.015 ±0.001 


10 


KERRDISK 


EW 


0.154 




0.153 


11 


IREFLECT 


Xpe 


5.34 ±0.15 


3.93 ±0.52 


4.72 ±0.12 


12 


IREFLECT 


S 


0.98 ±0.04 


1.00 ±0.10 


0.94 ±0.03 


13 


IREFLECT 




153.1 ± 15.7 


189.6 ±23.7 


153.3 ±11.6 



14 ^ 1.24(561/454) 0.82(297/361) 1.25(616/493) 



^ The parameter set is exactly the same as for Table 2. Case 1: Standard 0.7^5 keV fit results for our adopted model, which are copied from column SPl in 
Table 2. Case 2: Fit to SPl over the range 0.7^5 keV excluding 5.0-10.0 keV (i.e., excluding the Fe Ka line and Fe absorption edge). The value of q is that 
returned by KERRCONV. Case 3: Fit to SPl including HEXTE data that covers the full energy range from 0.7 to 150 keV. 
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Table 4 

The three nonrelativistic models 



Model 



NRl 
NR2 
NR3 



CRABCOR*CONST*TBABS(DISKBB+COMPBB+PEXRIV+GAUSSIAN) 
CRABCOR*CONST*TBABS(KDBLUR(g)(DISKBB+COMPBB+PEXRIV)+LAOR) 
CRABCOR*CONST*TBABS(KDBLUR0(SIMPLR0DISKBB+PEXRIV)+LAOR) 



APPENDIX 



A: PRELIMINARY NONRELATIVISTIC ANALYSIS 



With DISKBB as the principal component, we analyzed our three spectra in turn using three composite models of increasing 
sophistication, the Models NR1-NR3 Hsted in Table Al. In addition to diskbb, each model includes (i) a Compton power-law 
component, either COMPBB (Nishimura et al. 1986) or SIMPLR (Steiner et al. 2009b, 2011) convolved with DISKBB, and (ii) a 
reflection component PEXRiv (Magdziarz & Zdziarski 1995)^^, which models the absorption edges, plus an Fe Ka emission-line 
feature Gaussian or laor (Laor 1991). As in the case of our adopted model, also included are the three pre-factors crabcor, 
CONST and TBABS (see Section 3.2). In the case of Models NR2 and NR3, relativistic blurring effects are included using 
KDBLUR (Laor 1991). 

Model NRl: Cui et al. (1998) analyzed spectrum SPl using this model, and we do likewise. For the thermal component, 
our values of the parameters (disk temperature Tin and inner disk radius Ri^) are quite similar to those found by Cui et al.: 
respectively, Tin = (0.449 =b 0.005) keV (our Table 5) and Tin = (0.436 =b 0.004) keV (Table 2 in Cui et al. 1998); and 
R.^ = (2.54 zb 0.06)r^ versus (2.68 =b 0.06)r^ (for D = 1.86 kpc, M = 14.8 and i = 27.1 deg; Vg = GM/c^ = 21.9 km). 
For most of the other fit parameters, however, the differences between our results (Table A2) and theirs are quite significant: For 
example, for the width of the Fe Ka line we find cr = (1.22 ± 0.04) keV versus cr = (0.35 ± 0.07) keV, and for the ionization 
parameter ^ ~ versus ^ > 11911. We attribute these differences to our use of response files that have been significantly updated 
and improved (Section 2). These new response files not only allowed us to extend the upper energy bound to 45 keV (compared 
to the 30 keV bound used by Cui et al.), but they also greatly improved the fit: = 0-88 versus = 1.47. Furthermore, this 
comparison is understated because Cui et al. included 1% systematic error in the PCU count rates whereas we used 0.5%; for a 
systematic error of 1%, our = 0.68. 

Applying Model NRl to SP2 and SP3 also gives very good fits (Table A2, columns 5 and 6). However, comparing the results 
for all three spectra one finds a wide variation in the values of the disk temperature Tin (0.27-0.54 keV) and inner disk radius R-^^ 
(1.95-8.30 Vg). (We focus here on the parameters for the thermal component because this component ultimately delivers our key 
result, the spin of the black hole.) We find the performance of this model unsatisfactory for two principal reasons. First, it returns 
unlikely and near-zero values of the ionization parameter ^ (which we have here fixed to zero); this is probably because both the 
reflection component and line feature are improperly modeled (see below). Secondly, this model implies an unrealistically strong 
line feature (EW = 0.4-0.5 keV). Therefore, we now consider an improved model. 

Model NR2: For the line feature, we replace the symmetric Gaussian profile with a skewed, relativistic profile via the LAOR 
model, while relativistically blurring the three other additive model components using kdblur (Table Al). The LAOR model 
and its companion convolution model kdblur have the serious limitation of assuming a fixed and extreme value of the spin pa- 
rameter, a* = 0.998, and they therefore do not give a proper description of black holes with moderate or low spins. Nevertheless, 
we consider these models adequate for the present purposes because Cygnus X-1 is a rapidly spinning black hole and because the 
Fe Ka line is only a cosmetic feature relative to the thermal continuum (Section 4). The merit of using LAOR/kdblur here is that 
the model is simple and the values of the line strength and ionization parameter returned by the fits are reasonable (Table A3). 
In the end, however, as in the case of the previous model, we find Model NR2 wanting because the parameters it returns for 
the crucial thermal component (Tin and R{rd differ quite significantly for the three spectra. We now discuss our most physical 
nonrelativistic model. 

Model NR3: In addition to our criticisms of Models NRl and NR2, we find these models structurally unsatisfactory because 
they incorporate two different and discrepant thermal temperatures: DISKBB 's Tin ^ 0.5 keV and COMPTT's Tq ~ 1.0 keV. In 
Model NR3, we solve this problem by replacing COMPBB with SIMPLR, which generates the Compton power-law component via 
a convolution by operating on an arbitrary source of seed photons, which in this case is the thermal component DISKBB prior to 
Compton scattering. Meanwhile, we retain the line model LAOR and the convolution model kdblur. Model NR3 is specified 
in Table 4, and our fit results are given in Table 7. The fits are good and very comparable to those obtained using the other two 
models, even though Model NR3 uses one less fit parameter. We note that achieving these fits requires iron abundances that 
are ~ 3.2 — 5.5 times solar. The benefit of using this more self-consistent model is that it harmonizes the fitted values of the 
parameters of the thermal component: Tin = 0.532, 0.539 and 0.543 keV and i^in = 2.06, 2.30 and 2.01r^ for SPl, SP2 and 
SP3, respectively. Meanwhile, the smallness of the inner disk radius i^in is suggestive of the high spin revealed by our relativistic 
analysis (see Section 3). 



Four Preliminary Models, R1-R4: We now consider a progression of four models that are all built around our relativistic 
disk model KERRBB2 (Section 3.2). These models progress toward our adopted model (Section 3.2) in the sense that Model Rl 

As for our adopted model, we compute only the reflected component by Section 3 2) 
restricting the fit parameter s to negative values, which we here fix to -1 (see 



B: RELATIVISTIC ANALYSIS 



The Extreme Spin of Cygnus X-1 



17 



Table 5 

Fit results for Model NRl' 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 


DISKBB 


Tin 


0.449 ± 0.005 


0.537 ±0.010 


0.270 ±0.013 


2 


DISKBB 


^in 


2.54 ±0.06 


1.95 ±0.10 


8.30 ±0.98 


3 


CONST 


- 


1.012 ±0.002 


1.071 ±0.013 


1.037 ±0.004 


4 


TBABS 




0.675 ± 0.008 


0.697 ± 0.018 


0.862 ± 0.022 


5 


GAUSSIAN 


El 


6.40 ±0.14 


6.40 ±0.14 


6.40 ±0.12 


6 


GAUSSIAN 


a 


1.22 ±0.04 


1.02 ±0.07 


1.02 ±0.04 


7 


GAUSSIAN 


A^L 


0.047 ± 0.005 


0.078 ±0.011 


0.067 ± 0.006 


8 


GAUSSIAN 


EW 


0.488 


0.401 


0.455 


9 


COMPBB 


To 


0.777 ± 0.027 


0.947 ± 0.049 


0.432 ± 0.021 


10 


COMPBB 


Te 


29.00 ± 0.47 


21.33 ± 1.24 


23.15 ±0.64 


11 


COMPBB 


T 


1.08 ±0.02 


1.24 ±0.07 


0.80 ±0.05 


12 


COMPBB 


N 


4887.9 ± 741.6 


4074.1 ±978.0 


164572.0 ±48271.7 


13 


PEXRIV 




4.18 ±0.12 


2.08 ±0.25 


2.29 ±0.08 


14 


PEXRIV 


N 


9.07 ±0.35 


13.17 ±0.84 


14.44 ± 0.41 


15 






0.88(397/453) 


1.13(326/289) 


1.07(654/613) 



^ For the model components given, the parameters from top to bottom are: (1) inner disk temperature in keV; (2) inner disk radius in units of rg = GM/c^ = 
21.9 km for M = I4.8M0; (3) detector normalization constant relative to RXTE PCU2; (4) hydrogen column density in units of cm~^; (5) central line energy 
in keV; (6) line width in keV; (7) line flux in units of photons cm~^ s~^; (8) equivalent width of line in keV; (9) blackbody temperature in keV; (10) electron 
temperature of corona in keV; (11) optical depth of corona; (12) normalization constant; (13) iron abundance relative to solar; (14) normalization constant; (15) 
Reduced chi-square, total chi-square and degrees of freedom, respectively. Details for the reflection component PEXRiv: Apart from iron, the metal abundances 
are solar; the photon index is fixed at F = 2.5, the value determined for Model NR3; the reflection scaling factor s is fixed to -1; and the ionization parameter is 
fixed, ^ = (whereas, if fitted, ^ 0). 
^ The line central energy is pegged in the fit. 



Table 6 

Fit results for Model NR2^ 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 


DISKBB 


Tin 


0.509 ± 0.003 


0.616 ±0.011 


0.446 ± 0.002 


2 


DISKBB 


^in 


2.09 ± 0.02 


1.61 ±0.05 


2.40 ±0.01 


3 


CONST 




1.014 ±0.002 


1.083 ±0.011 


1.034 ±0.004 


4 


TBABS 


Nn 


0.701 ±0.007 


0.801 ±0.016 


0.717 ±0.006 


5 


LAOR 


El. 


6.51 ±0.02 


6.49 ± 0.04 


6.49 ± 0.02 


6 


LAOR 


Q 


2.37 ±0.07 


2.54 ±0.11 


2.39 ± 0.04 


7 


LAOR 


A^L 


0.012 ±0.001 


0.035 ± 0.003 


0.032 ± 0.001 


8 


LAOR 


EW 


0.125 


0.176 


0.219 


9 


COMPBB 


To 


1.060 ±0.015 


1.279 ±0.065 


0.522 ±0.001 


10 


COMPBB 


Te 


33.64 ±0.64 


23.23 ±1.56 


27.20 ± 0.42 


11 


COMPBB 


T 


1.12 ±0.02 


1.22 ±0.08 


0.87 ±0.02 


12 


COMPBB 


N 


1404.2 ±83.0 


1165.7 ±270.6 


59084.1 ±426.3 


13 


PEXRIV 




2.96 ±0.37 


2.31 ±0.38 


3.63 ±0.17 


14 


PEXRIV 




1168.0 ±115.4 


748.5 ± 177.2 


1647.7 ± 115.8 


15 


PEXRIV 


N 


4.91 ±0.32 


14.51 ±0.98 


10.33 ± 0.35 


16 




xi 


1.05(474/452) 


1.16(333/288) 


1.13(692/612) 



^ Layout and parameter definitions very similar to Table 3, with two exceptions: (1) LAOR model substituted for GAUSSIAN and, correspondingly, the emissivity 
index q is given in place of the line width cr; and (2) the ionization parameter ^ is now a fit parameter rather than fixed to zero. 



is the most primitive and our adopted model is the most physically realistic. All four of these models and our adopted model give 
very similar results for the parameter of interest, namely a*, indicating that our key result, the extreme spin of Cygnus X-1, is 
quite insensitive to the details of the analysis. 

The four preliminary models R1-R4 are defined in Table 8. Every individual component in all four models has already been 
described either in Section 3.2 or in Appendix A. There are only two combinations of model components that are new^ and that 
have not been used elsewhere, namely, the reflection components in Models R3 and R4. Their core power-law components are 
respectively COMPBB and SIMPLC, each of which is convolved in turn with IREFLECT and KDBLUR. 

The results for Models R1-R4 are presented respectively in Tables 9-12. Comparing the results for Models Rl and R2, 
which both employ the widely-used reflection model PEXRIV, Model R2 is preferred because it gives very comparable values of 
reduced chi-square using two fewer fit parameters, and its reflection component SIMPLR self-consistently generates the power- 
law with no need for the second thermal component required by COMPBB. Furthermore, in the case of the nonrelativistic analysis 
(Appendix A), it was shown that simplr delivers consistent values of Ri^ and Tin, which COMPBB fails to do. 

Tables 1 1 and 12 summarize our results for Models R3 and R4, which use the convolution reflection model ireflect in place 
of the additive model PEXRiv. While both reflection models are based on the same physics, the virtue of ireflect is that it can 
take any shape of spectrum as input, whereas PEXRIV requires a power-law input spectrum. In considering Models R1-R4, we 
find Model R4 to be superior because it fits the data as well as the other models while using the fewest parameters, and it is the 
most self-consistent (see Section 3.2). 

In a final step, as described in Section 3.2, we arrive at our adopted model by replacing the flawed model components LAOR and 
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Table 7 

Fit results for Model NR3^ 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 


DISKBB 


Tin 


0.532 ± 0.001 


0.539 ± 0.002 


0.543 ±0.001 


2 


DISKBB 


^in 


2.06 ±0.01 


2.30 ±0.01 


2.01 ±0.01 


3 


CONST 




1.018 ± 0.002 


1.078 ± 0.010 


1.0295 ± 0.0038 


4 


TBABS 


Nn 


0.680 ± 0.004 


0.838 ± 0.009 


0.699 ±0.004 


5 


SIMPLR 


r 


2.206 ±0.011 


2.503 ±0.013 


2.486 ±0.005 


6 


SIMPLR 


/sc 


0.176 ±0.002 


0.344 ± 0.007 


0.317 ±0.003 


7 


LAOR 


El 


6.56 ± 0.02 


6.44 ± 0.04 


6.49 ± 0.03 


8 


LAOR 


q 


2.88 ±0.01 


2.36 ± 0.06 


2.38 ±0.03 


9 


LAOR 


A^L 


0.018 ±0.001 


0.033 ± 0.002 


0.030 ± 0.001 


10 


LAOR 


EW 


0.189 


0.164 


0.208 


11 


PEXRIV 




5.49 ±0.17 


3.24 ±0.13 


3.22 ±0.14 


12 


PEXRIV 




1147.3 ± 100.4 


906.5 ±81.5 


1323.5 ±0.0 


13 


PEXRIV 


N 


3.78 ±0.29 


15.38 ±1.17 


9.52 ±0.30 


14 




xi 


1.26(570/454) 


1.27(367/290) 


1.12(685/614) 



^ For definitions of most of the parameters, see Tables 5 and 6. The model components and parameters are the same as for Table 6 with one difference: The 
convolution model SIMPLR with two parameters has been substituted for the additive model COMPBB with four. The two parameters of SIMPLR are the photon 
power-law index F and a normalization constant, which is the scattering fraction /sc- 



Table 8 

The four preliminary relativistic models 



Model 



Rl CRABCOR*CONST*TBABS(KERRBB2+LAOR+KDBLUR0(PEXRIV+COMPBB)) 

R2 CRABCOR*CONST*TBABS(SIMPLR(8)KERRBB2+LAOR+KDBLUR(g)PEXRIV) 

R3 CRABCOR*CONST*TBABS(KERRBB2+LAOR+KDBLUR(g)(lREFLECT0COMPBB)) 

R4 CRABCOR*CONST*TBABS(SIMPLR(g)KERRBB2+LAOR+KDBLUR(g)(lREFLECT(g)SIMPLC)) 



Table 9 

Fit results for Model Rl' 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 


KERRBB2 
KERRBB2 

const 

TBABS 
LAOR 
LAOR 
LAOR 
LAOR 
COMPBB 
COMPBB 
COMPBB 
COMPBB 
PEXRIV 
PEXRIV 
PEXRIV 


a* 
M 

TVh 

^L 

Q 

Ni. 
EW 
To 

Te 
T 

N 
N 


0.102 ±0.007 

1.000 ±0.002 

0.760 ± 0.008 
6.45 ± 0.02 
2.07 ±0.10 

0.011 ±0.001 
0.106 

1.023 ± 0.005 
24.23 ± 0.47 
1.11 ±0.01 

1582.7 ±24.5 
3.01 ±0.30 

751.3 ± 104.9 
4.87 ±0.31 


U.yyyy_Q QQQ3 

0.143 ±0.030 
1.004 ±0.014 
0.893 ± 0.022 
6.42 ± 0.04 
2.02 ±0.20 
0.026 ± 0.003 

0.131 
0.975 ± 0.028 
21.41 ± 1.06 
1.32 ±0.07 
3851.8 ± 504.1 

2.87 ±0.38 
652.1 ± 169.2 
12.03 ± 1.07 


U.yyOy_Q QQ21 

0.094 ± 0.010 
1.016 ±0.005 
0.759 ±0.010 
6.48 ± 0.03 
2.33 ±0.05 
0.030 ± 0.001 

0.208 
0.559 ± 0.002 
24.23 ±0.57 
1.05 ±0.03 
39574.0 ± 515.0 
3.63 ±0.18 
1226.3 ±0.1 
10.87 ±0.40 
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1.04(470/452) 


1.26(364/288) 


1.16(708/612) 


^ For the model components given, the parameters from top to bottom are: 


(1) spin parameter; (2) mass accretion rate in units of 10^^ 


g s~^; (3) detector 



normalization constant relative to RXTE PCU2; (4) hydrogen column density in units of cm ; (5) central line energy in keV; (6) emissivity index; (7) line 
flux in units of photons cm~^ s~^; (8) equivalent width of line in keV; (9) blackbody temperature in keV; (10) electron temperature of corona in keV; (11) 
optical depth of corona; (12) normalization constant; (13) iron abundance relative to solar; (14) ionization parameter; (15) normalization constant; (16) Reduced 
chi- square, total chi-square and degrees of freedom, respectively. Details for the reflection component PEXRIV: Apart from iron, the metal abundances are solar, 
and the photon index is fixed at F = 2.5, which is the value determined for our adopted model (Table 2). 

KDBLUR (see Appendix A) by kerrdisk and kerrconv. The virtue of kerrdisk/kerrconv compared to laor/kdblur 
is that the spin parameter is free, allowing the GR metric to be calculated properly. We use a single, linked value of a* in fitting 
kerrdisk/kerrconv and kerrbb2. An inspection of Tables B2-B5 and Table 3 shows that the key quantity, namely the spin 
parameter, is precisely determined and near unity for Models R1-R4 and our adopted model. This indicates that the details of 
how one models the power-law and reflected components has little affect on the determination of a*. 
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Table 10 

Fit results for Model R2' 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 


KERRBB2 
KERRBB2 
const 
TBABS 
SIMPLR 
SIMPLR 
LAOR 
LAOR 
LAOR 
LAOR 
PEXRIV 
PEXRIV 
PEXRIV 


a* 
M 

Nn 

r 

fsc 

^L 

q 

A^L 

EW 
N 


9989+°-°°°^ 

U.»»Oy_Q QQQ3 

0.116 ±0.003 
1.006 ± 0.002 
0.758 ± 0.006 
2.295 ±0.011 
0.179 ±0.002 
6.54 ± 0.03 
2.87 ±0.02 
0.014 ±0.001 
0.148 
6.23 ±0.17 
726.4 ±39.5 
6.91 ±0.44 


0.99981q QQ3g 

0.146 ±0.025 
1.042 ± 0.013 
0.841 ±0.018 
2.507 ±0.012 
0.328 ± 0.009 

6.41 ± 0.04 

2.27 ±0.07 
0.032 ± 0.002 
0.154 

3.06 ±0.12 
981.7 ± 103.4 

15.23 ±1.01 


QQQQ+0 

U.»yyy_Q QQQJ^ 

0.124 ±0.010 
0.972 ± 0.004 
0.764 ±0.010 
2.525 ± 0.008 
0.299 ± 0.003 

6.51 ±0.03 

2.41 ±0.05 
0.028 ± 0.001 
0.187 

3.60 ±0.15 
1163.5 ±73.4 

12.28 ±0.53 


14 




xi 


1.24(561/454) 


1.28(371/290) 


1.18(725/614) 



^ For definitions of most of the parameters, see Table 9. Two distinctions between this table and Table 9: (1) The convolution model SIMPLR with two parameters 
has been substituted for the additive model compbb with four. The two parameters of SIMPLR are the photon power-law index F and a normalization constant, 
which is the scattering fraction /sc- (2) The photon index in PEXRIV is not fixed, rather it is linked to the photon index in SIMPLR. 



Table 11 

Fit results for Model R3^ 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 


KERRBB2 
KERRBB2 

const 

TBABS 
LAOR 
LAOR 
LAOR 
LAOR 
COMPBB 
COMPBB 
COMPBB 
COMPBB 
IREFLECT 
IREFLECT 
IREFLECT 


a* 
M 

^L 

q 

A^L 

EW 
To 

Te 
T 

N 

Xpe 
S 

i 


QQ8Q+0-0003 

0.110 ±0.004 
1.001 ±0.002 
0.733 ± 0.006 
6.43 ±0.02 
1.78 ±0.14 
0.011 ±0.001 

0.109 
0.942 ± 0.010 
26.25 ± 0.43 

1.19 ±0.02 
2212.8 ±99.1 

2.20 ±0.19 
0.36 ± 0.02 
84.2 ± 9.0 


QQQQ+0 

U.»»»»_Q QQQ3 

0.145 ±0.031 
0.984 ±0.015 
0.820 ± 0.025 
6.42 ± 0.04 
1.83 ±0.28 
0.024 ± 0.003 

0.117 
1.004 ±0.017 
22.85 ±0.61 
1.24 ±0.04 
3301.5 ±234.5 
2.02 ± 0.26 
0.72 ± 0.03 
59.2 ± 12.5 


0.098 ±0.008 
1.026 ±0.004 
0.670 ±0.011 
6.45 ± 0.03 
2.26 ±0.06 
0.030 ± 0.001 

0.204 
0.518 ±0.001 
26.25 ± 0.48 
0.90 ± 0.03 
59300.0 ± 532.3 
3.51 ±0.18 
0.72 ± 0.03 
65.5 ±6.8 


16 




Xv 


1.07(483/452) 


1.25(358/288) 


1.14(695/612) 



^ Here, we have substituted the convolution model IREFLECT for the additive refiection model PEXRiv. These models have two parameters in common, the iron 
abundance and ionization parameter, and one that differs, namely, the normalization of PEXRIV is replaced by the reflection scaling factor s. For the definitions 
of all other parameters, see Tables 9 and 10. 

^ In applying IREFLECT, the reflection scafing factor s is negative; here we give absolute values of s. 



Table 12 

Fit results for Model R4^ 



Number 


Model 


Parameter 


SPl 


SP2 


SP3 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 


KERRBB2 
KERRBB2 

const 

TBABS 
SIMPLR 
SIMPLR 
LAOR 
LAOR 
LAOR 
LAOR 
IREFLECT 
IREFLECT 
IREFLECT 


a* 
M 

r 

fsc 
El. 

q 

Ni. 
EW 

Xpe 
S 


0.9987i°:ro^ 
0.115 ±0.004 
1.000 ±0.002 
0.704 ±0.006 
2.284 ± 0.010 
0.227 ± 0.006 
6.54 ± 0.03 
2.85 ± 0.02 
0.016 ± 0.001 
0.161 
5.33 ±0.15 
1.00 ±0.04 
148.0 ± 14.7 


QQQ7+0.0001 

u.yyy/_Q QQ2g 
0.143 ±0.026 
1.037 ±0.013 
0.761 ± 0.022 
2.525 ±0.012 
0.329 ±0.013 
6.42 ± 0.04 
2.26 ±0.07 
0.030 ± 0.002 
0.150 
3.16 ±0.12 
0.89 ±0.03 
53.6 ±7.9 


0.9999ig.oooi 
0.122 ±0.011 
0.726 ± 0.003 
0.687 ±0.010 
2.551 ±0.011 
0.308 ± 0.006 
6.48 ± 0.04 
2.33 ± 0.07 
0.025 ± 0.002 
0.165 
3.63 ±0.14 
0.95 ± 0.03 
60.8 ±8.2 


14 




Xv 


1.20(545/454) 


1.28(370/290) 


1.17(720/614) 



^ Same model components and parameters as in Table 11, except that we replace the power-law model compbb with the convolution model simper, thereby 
reducing the four parameters of the former model component to two, namely the photon power-law index F and the normalization constant /sc- 



